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of the transformation 


This description 


telephone industry in the electronic 


based upon the development in the bries 
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however, that the experience 


Co.is of general application. Probably 


noanethel its transformation has been more 


e 
drastic than elsewhere since LME hardly had 
any markets protected from 


with old 


competition where It 


could continue the technology. 


One may distinguish three functions 


telephone network: the conversion function, 


the transmission function, and the switching 


function. The telephone industry designs and 


manufactures equipment for these functions, 
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the conversion function — telephon sets 


In telephone sets a conversion of acoustic to 


electric energy takes place, and vice versa. 


The telephone transmits and receives speech 
currents. 


transmission function — cable and 


By means of cable, wire, carrier systems, the 


electric speech currents are transmitted from 


one place to another. 


the switching function — telephone 
exchanges 
subscribers 


The telephone exchanges con the 


together according to 
A telephone exchange alao 
the trè 


concentrating 


lisation of 
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With regard to Fig. 1 and 2 "Key achieve- 


ments'' some comments are made, 


Telephone sets have remained remarkably 


unchanged for nearly a century. Their components 


EL 
Th 
JA 
[D 
[u 
Ls | 
hl] 


result of new produc- 


tion methods and new materials. In the sixties 
and seventies, however, a revolution oecurred. 


underwent radical changes. 


ld carbon granule microphone has been 


better sound repro= 


and life properties. The technigue now 


ectromagnetie, piezoelectrical 


microphones ete. On the other hand 


do not have the fine property of the 
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carbon granule the inherent gain. 


mierophone, 
The new microphones require external gain,which 


effected with semiconductor circuits. 


dial is being more and more superseded by 


In some cases new acoustic 


signal 
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very 
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advances have been most striking. The analog 


systems predominated until quite recently; the 


investments in digital equipment. 


main 


Over the years some different epochs are 


suishable as regard 


s components 
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KEY ACHIEVEMENTS KEY ACHIEVEMENTS 


1844 The telegraph (Morse) 

1865 ITU 

1876 The telephone (Bell) 

1889 The Strowger selector 

1901 Transatlantic wireless telegraphy (Marconi) 
1906 The triod (Lee de Forest) 

1910 Teletraffic theory (Erlang) 


1923 Telephotography (lves) 
The iconoscope (Zworykin) 


1927 The feed back amplifier (Black) 


1930 The coaxial cable 
Telex (Siemens) 


1937 Pulse code modulation (Reeves) 


1942 Information theory (Shannon) 

1946 ENIAC (Maucly & Eckert) 

1947 The transistor (Bardeen, Brattain & Shockley) 
1958 The integrated circuit (T | & Fairchild) 

1962 Telstar (ATT) 


1965 Intelsat in syncronous orbit (Rosen). 
Stored program control ESS 1 (ATT) 


1966 Conception of step index optical fibre (Kao) 
1970 Laboratory tests with first practical optical fibre (Maurer) 
1971 The microprocessor (Intel) 


1976 Fiber optical transmission in operation (ATT) 
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Fig. History of Telecommunication — key Fig.2 Hitstory of Telecommunication — key 
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Optical fibres becoming increasingly common and will be the 


Switching technology has passed through three predominant practice 


in the future. The other 
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This is particularly the case with L.M.Ericsson, Fig.3 shows the development so far and what 
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equipment both from a technical point of view another way 
and for produetion volume. Â erossbar exchange contains per subscriber 


Before going in detail into the changes line 4-10 relays and out 1.4 crossbar verticals 
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in our structure a few worde will be mentioned Â electronic exchange contains more than 10 000 
nr 


about what the new electronics involves. transistors per line! 


A characteristic of the new electronica is A part of the difference is due to the fact that 
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7 . . Farscas ts Versus Reality ‘Ericsson s shift from 
Number of circuit elements per component te. mechanical to Electronic Telephone Exchanges 
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Fig.d Shift from Eleectro-mechanical to Electro- 


nie Telephone Exchanges — Forecasts 


versus reality 


How, then, has the conversion to electronics at, at most, the same price as they can be pur- 
worked out in practice? Fig.4 shows that the chased. To start with, therefore, two categories 


production of electromechanical systems reached come under consideration for manufacture within 
its peak at the end of 1975. It then fell very the Group: a) components produced in large 
sharply, Tw (1980) being stabilized at about volumes, and 5) special purpose components, i.e. 
one quarter of the peak level, Electronic produc- components designed exclusively for the 

tion increased at the same time, but, to start telephone industry. 

with, much more slowly. Now, however, it would This extensive reorganigzation of the produc- 
seem to have reached roughly the same value as tion brings with it many problems. These may be 


that of electromechanical production in 
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said to exist three spheres: production 


The picture is blurred by the fact total engineering, training of staff, and psychological 
sales fell in the years 1976-f/. Of greater problems, The production technigue naturally 
interest is to see how the forecasts changed requires very thorough preparation. A careful 
successively towards a higher proportion oi study must be made of the processes one wishes 
electronic and lower proportion of eleectromecha- to adapt. 
nical production. This, of course, impose very It is important to acquire the correct maehines 
exacting requirements on flexibility in the for the anticipated production volumes =— compar- 
production, atively large resources must be allocated to 
Fig.5 shows the composition of the produc- testing. These problems are well known in all 
tion costa for different kinds of materials. On new production and are solved by input of 
ion to AXE analog there is a reduction engineering work and capital. 
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In conjunction with the present shift to digital would be more troublesome, Shop-floor workers 


AXE there is a further reduction to 51%. At the must ask themselves: is the company management 
same time Ericsson becomes a large purchaser of right in going in for this technique? 
electronic circuits on the world market. As & Can I hold dewn a job in electronics production? 


result of the drop in value-added in Ericsson's Shall I be discharged owing to lack of work,shall 
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HOW ELECTRONICS HAS CHANGED TiE STRUCTURE 
OF MANUFACTURING COSTS - 
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Fig.5 Change of structure of manufacturing 


costs 


It turned out to be fairly easy to convince 
workers of the necessity to go for electronics. 
Swedish workers have traditionally had a positive 
attitude to technical innovations even if their 
jobs appear to be treatened by them. So, too, 
in this case: they even accept a reduced level 
of earnings (for workers who, through skill and 
hard work, had been able to earn good piecework 
wages). In some cases, perhaps, it was more 
difficult to convince white-collar workers of 
the necessity for changing to electronics. 

Some of them had become specialists in particu- 
lar fields and saw that their expert knowledge 
would be of less value. 

Many workers (of both sexes) who had reach 
ed a certain age were worried as to whether they 
could successfully tackle the new type of work. 
There were even people ín supervisory positions 
who doubted their ability. But the fears proved 
to be exaggerated. After a fairly short period 
of adjustment the workers have proved very 
capable in the new production. The impairment 
of eyresight and dexterity which we all suffer 
to different degrees after 50 years of age has 
not been any hindrance. 

There was, of course, an alternative to 
retraining of the workers for electronics produc- 
tion and that was to take on and train new 
people. But this would have meant discharging a 
large number of the earlier workers — a course 
that was impossible for an employer with social 
responsibility, and also very difficult on legal 
grounds. Instead of building new factories for 
electronics with an entirely new work force, 


it was decided to convert some of the old 
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_15,000 


10,000 


White collar 
5,000 (incl. technical staff) 
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company's labour forca in Swenden | axchoting 
op layars on military contracts, and adjusted to 
allem lor partial corporate restructuring during tha period 
Fig.6 Change of manpower mix over the years 


1975-78 


factories, wholly or partly, and to retain the 
existing work force. In this way a conflict with 
the trade unions could be avoided and the good 
spirit in the company could be maintained. 

On the engineering side we had some diffi- 
eulty in recruiting software people. There is a 
shortage of such people throughout the world. 
Ericsson decided to decentralize this work 
wherever possible. Thus there are groups of 
people in many countries engaged on software 
work for delivery of equipment to their own 
countries. This is advantageous, as the work is 
usually done in cooperation with the customer, 
the local telephone administration. It also 
brings about a spread of competeney within 
the Group, to which we give high priority. 

The changes in switching technology are 
those that have meant most for the structure of 
the telephone industry. 

Innovations on the telephone set side, involving 
a consderably larger proportion of electronics, 
have a similar trend; but as the volume is 
comparatively small, the consequences are not 
so radical. Telephone sets contain both analog 
and digital technigue. We are probably quite a 
long way off from a general digitalization of 
telephone sets, Developments in office organi 
zation and communication equipment for office 
may, however, place new requiremente on 
electronic production capacity. 

Transmission technology has always been 
electronic apart from the physical circuits.Since 


the latter will in due course become electronic to 
come extent, there will be — here too — an 


increased emphasis on electronics. 


PUBLIC NETWORK INVESTMENTS 
TELEPHONY 
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Fig. Public network investments in telephony 


\not commented upon in paper) 


Finally a few worda about the advantages and savings can be made also since electronic | 
disadvantages cf the electronic exchanges from exchanges occupy less space. On the other hand 
the user's standpoint. they require more energy for their operation. 
As regards those whom the exchanges serve, the From the maintenance point of view they are 
subscribers, the advantages may be said to be rather better, but more important perhaps is 
fairly limited. Telephoning takes place roughly tnat administrations find it easier to recruit 
as before. Some new facilities can be introduced, statt for eleetronies than tor electromechanics. 
such aa abbreviated dialling (i.e. each sub- Chief ly, however, the administrations place a 
scriber can call his commonest "correspondents'' value on the greater flexibility of the 
by dailling three-digit numbers), automatic electronie exchanges. They expect, therefore, 
transfer, automatic alarm clock, ete., but that adaptation to new requirements will be 
experience hitherto shows their use to be limi- easier. The way is also opened to 100% disital 
ted. telephone plant with integrated switching and 
With 100% digitalization the transmission charac- transmission. 
teristics can be standardigzed within narrow All of this gives administrations the hope of 
limite, but the subscribers will notice this some reduction of costs, at least in the long 
only in special cases. For the adiminiatrations term, which will also benefit subscribers. 
the situation is rather different. The electronic Even if the local network can also be made 
exchanges are today hardly cheaper than the cheaper, since the exchanges can be divided up, 
electromechanical, but there are good prospects so reducing the average length of subscriber 


that their prices will fall fail in the future. lines, the gains will be moderate. 


Voordracht gehouden op 11 september 1980 in de Afdeling 
Elektrotechniek THD tijdens een gemeenschappelijke ver- 
adering van het NERG (nr.291), de Sectie voor Telecom=- 


municatietechniek KIvIl, en de Benelux sectie IEEE. 


NEDERLANDS ELEKTRONICA- EN RADIOGENOOTSCHAP 
Regionale Bijeenkomst Rotterdam 
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In het kader van het 60-jarig bestaan van het NERG worden er enkele korte regionale bijeenkomsten 
georganiseerd, waarin wordt stilgestaan bij de voortgaande toepassing van de elektronica in andere 
vakgebieden en de enorme maatschappelijke betekenis daarvan. Het programma van een regionale 
bijeenkomst is geheel anders van opzet dan dat van de gangbare werkvergaderingen; niet alleen 
NERG-leden, maar ook hun familieleden en andere geïnteresseerden krijgen in kort bestek een 
overzicht van bepaalde ontwikkelingen. Verder zal er een ongedwongen ontmoeting kunnen zijn met 
vertegenwoordigers van andere vakgebieden, die op grote schaal gebruik maken van de mogelijk- 
heden van de moderne elektronica. 

Op 8 november 1895 ontdekte Wilhelm Conrad Röntgen, die net buiten onze landsgrenzen geboren is 
maar in Nederland opgegroeid, een nieuwe stralensoort. Wat er uit deze ontdekking nog eens zou 
voortvloeien kon hij niet overzien, maar ruimschoots genoeg voor een feestelijke bijeenkomst 
met als thema: 


85 jaar 
Röntgenstraling 


60 jaar 
Elektronica 


van ontdekking tot computer-tomogratie 


Deze bijeenkomst wordt gehouden in de conferentiezaal van het Academisch Ziekenhuis „Dijkzigt” 
te Rotterdam op zaterdagmorgen 8 november 1980 van 10.00 tot 12.30 uur. Gastheren zijn 
Prof. K. Hoornstra (Erasmus Universiteit) en Ir. F. v. d. Meer (Dijkzigt-Ziekenhuis). 

Er zijn belangrijke ontwikkelingen gaande op het gebied van de computer-tomografie (uitgaande 
van röntgenstraling of kernspinresonantie) binnen Philips Medical Systems; over deze ontwikke- 
lingen zal Ir. G. J. A. Arink wat vertellen. Voor deze bijeenkomst is de keus gevallen op een 
zaterdag omdat het op die dag alleen mogelijk is de apparatuur te bezichtigen. 

Parkeerruimte is op die dag doorgaans ook geen probleem. Het Academisch Ziekenhuis Dijkzigt is 
vanaf het N.S. Station Rotterdam C.S. snel te bereiken met tramlijn 4. 

Gezien de beperkte plaatsruimte en de te verwachten grote belangstelling wordt tijdige aanmelding 
verzocht door middel van de aangehechte kaart, die duidelijk moet worden ingevuld en gefrankeerd 
met een postzegel van 45 cent. Wie na aanmelding verhinderd mocht zijn dient dit even door te 
geven aan Ir. F. v. d. Meer, tel. 010-6535987. 


Namens het voorbereidingscomité, 
Rotterdam, oktober 1980. e W. Herstel. 


KOMPONENTEN FÜR DIGITALE VIDEOSIGNALVERARBEITUNG IN NMOS-TECHNOLOGIE 


P. Draheim 
VALVO — RHW — Philips GmbH Hamburg 


Der Bericht behandelt die Entwicklung von Kom- 
ponenten für wvideofähige A/D- und D/A-Wandler 
in NMOS-Technik. Es wird zunächst das unter 
anderem angestrebte Ziel eines 8 Bit A/D-Wand- 
lers nach dem Pipelining-Verfahren dargestellt 
und danach die schon realisierten Komponenten 
sowie ein 8 Bic D/A-Wandler. 

In Abb. 1 ist das Bloekschaltbild des 
8 Bit-A/D-Wandlers nach dem Pipelining-Verfah- 
ren zu senen. Das analoge Eingangssignal wird 
Einer Track + Hold-Schaltung zugeführt und von 
dort einem d-Bit Grobkonverter mit 15 Kompara= 
toren. Hinter dem Grobkonverter liegt ein Sig- 
nal im l-aus-16 Code vor, das in einen Binär- 
Code umgewandelt wird. Die in dem Grobkonver=- 
ter erzeugten, binär codierten 4 Most-Signifi- 
cant-Bit (MSB's, B4-B7) werden danach in einem 


Latch zwischengespeichert. 
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Abb. 1: 8-Bit A/D-Konverter nacn dem Pipe- 


lining Verfahren 


Weiterhin wird das l-aus-16 codierte Sig- 
nal hinter dem Grobkonverter einem D/A-Wandler 
zugeführt und das rückgewandelte Signal gleich=- 
zeitig mit dem Eingangssignal aus der Track + 
Hold-Schaltunag in einem l6-fach Differenzver- 
stärker verstärkt auf den Eingang des 4-Bit- 
Feinkonverters gegeben. Der Feinkonverter ist 
identisch wie der Grobkonverter aufgebaut und 
er erzeugt die Least-Significant-Bits (LSB's, 
BO=B3} . 
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Ein solcher A/D-Wandler ist vom Bauelementeauf- 
wand gegenüber einem voll parallelen 8-Bit A/D- 
Wandler mit 255 Komparatoren sehr ökonomisch 
und hat die gleiche Wandlungsrate wie ein voll 
paralleler. Jedoch ist diese Entwicklung aus 
zwei Gesichtspunkten nicht unproblematisch. 

Erstens muA der 4 Bit-Grobkonverter eine 
Linearität von 8 Bit haben, da die in ihm er- 
zeugten Linearitätsabweichungen um den Faktor 
16 durch den Differenzverstärker verstärkt dem 
Feinkonverter zugeführt werden und dort zu ei- 
nem fehlerhaften Ergebnis führen. 

Zweitens darf der Differenzverstärker kei- 
nen Offset haben und die absolute Verstärkung 
muB auf 16 + l d‚h. + 6 % genau sein. 

Für die Fehlerkorrektur des Grobkonverters 
sind Fehlerkorrekturschaltungen entwickelt wor= 
den, die bei einer A/D-Wandler-Eingangsspannung 
von 3 Vsg, Komparatorfehler bis zu + 80 mV im 
Grobkonverter ausregeln können und somit eine 
ausreichende Linearität gewährleisten. 

Zur Reduktion des Offsets im Differenzver- 
stärker wurde ein getakteter Verstärker entwik= 
kelt und durch die Anwendung eines Doppeldiffe- 
renzverstärkersystems, in dem auch die Referenz- 
spannung für den Feinkonverter von der 
Verstärkung abhängt, die Abhängigkeit 
von der absoluten Verstärkung von + 6 % auf die 
relative Verstärkungsdifferenz von + 6 % zweier 
benachbarter Differenzverstärker reduziert. 
Dies ist in integrierter Technik gut darstell- 
bar. Die in dem Pipelining-System angestrebte 


Spezifikation sieht folgendermaben aus: 


_ Auf lösung 8 Bit monoton 

_ Taktrate 20 MHz 

_ Bandbreite 6 MHz (3 dB) 

je Eingangssignal 3 Vas 

an Versorgungsspannung + 12 V + 10 & 
+ SV tE 20 8 

_ Verlustleistung < 1 W 

- Ausgangs Level TTL 

= Takt Level TTL 


Wie man der Spezifikation entnehmen kann zielt 
diese Entwicklung auf eine Anwendung in der di- 
gitalen Signalverarbeitung von Fernsehempfangs- 


geräten. 


Als erste Basiszelle wurde ein 4 Bit A/D-Wand- 
ler mit 15 Komparatoren und einem eingebauten 
D/A-Wandler in einem NMOS-ProzeB mit kleinsten 
strukturen von 4 um entwickelt. Abb. 2 zeigt 
das Chipfoto. 


Die aktive Chipfläche (ohne Bond-Pads) 
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Abb. 2: Chipfoto des 4 


Aufgrund der guten Linearität (> 5 Bit) wurde 
im gleichen ProzeB mit den gleichen Komparato- 
ren ein 6 Bit paralleler A/D-Wandler mit 63 
Komparatoren entworfen. 

Abb. j zeigt, das Logikdiagramm mit den 
63 Komparatoren Kl-K63, den nachgeschalteten 
Flip-Flops, hinter denen ein n-aus-64 Code 
vorliegt, den NOR-Gattern zur Umcodterung in 
den l-aus-64 Code und dem ROM 63x6 zur Binär- 
Codierung.In Abb. 4 ist das Chipfoto zu sehen. 
Dieser 6 Bit-A/D-Wandler hatte eine aktive 


= 


ie Ë 
Fläche von 2 mm und die gemessenen Werte wa- 


ren: 

ie Auf lösung 6 Bit 

r Linearität > 6 Bit 
_ Taktrate >20 MHz 
5 Bandbreite 5 MHz 
ei Verlustleistung Os MW 


Die Linearität war ausreichend um zwei 6 Bit= 


A/D-Wandler parallel zu schalten und so einen 


war ca. 1 mm? und dieser A/D-Wandler hatte fol- 


gende Eigenschaften: 


mi Auf lösung A Bit 

a Linearität + Bit 

= Taktrate „20 MHz 

a Bandbreite 6 MHz (3 dB) 


ns Verlustleistung < O.l W 


„ner TE 


none 


Bit A/D-Konverters 


{ Bit parallelen A/D-Wandler auf zubauen. 

Neben den hier aufgeführten A/D-Wandler 
Entwicklungen wurde auch ein 8 Bit D/A-Wandler 
entworfen. Abb. 5 zeigt das Logikdiagramm. Das 
im Binär-Code vorliegende Eingangssignal wird 
aut 255 NOR-Gatter gegeben, die dieses Signal 
in ein l-aus-256 Code-Signal umcodiert. Dieses 
l-aus-256 Code-Signal wird dann getaktet 255 
“Transmission-Gates" zugeführt, die dann, je 
nachdeft welches Transmission-Gate von dem 1- 
aus=256 Signal angesprochen wird, eine Span- 
nung von dem Referenzspannungsteiler auf die 
Ausgangsleitung schalten. Der Referenzspan- 
nungsteiler besteht aus einer Poly-Silizium- 
Bahn mit 255 Abgriffen. Der absolute Wider- 
standswert dieses Spannungsteilers ist nicht 
von Bedeutung, sondern nur die relativen Wi-= 
derstandsverhältnisse zwischen den einzelnen 
Abariffen. 


In Abb. 6 ist das Chip-Foto dieses D/A- 


Wandlers zu sehen, der ebenfalls eine aktive 


5 
Chipfläche von 2 mm hatte. 


_é-bit-A/D-Wandler 


wehiehikjt. 
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Abb. 3: Logikdiagramm 
Die Linearität war besser als 1/2 LSB, die Abb. 4: Chipfoto eines 6 Bit A/D-Konverters 
Taktrate gröBber 20 MHz und die Verlustleistung 


kleiner O,3 Watt. 
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RELAXATIE-METHODEN IN DIGITALE BEELDSEGMENTATIE 


Ir. J.J. Gerbrands 


Technische Hogeschool Delft 


Afdeling der Elektrotechniek 


Vakgroep Informatietheorie 


Relaxation methods in digital image se 


entation. Probabilistic relaxation label procedures are discussed 


in the context of low-level image segmentation. Two examples are discussed in detail. The first is the re- 


inforcement of edges as proposed by Schachter et al. The second example deals with image segmentation by 


clustering in the feature domain followed by relaxation in the spatial domain. The features used are gray- 


level and texture. 


INLEIDING 


Tot het gebied van de digitale beeldverwerking kan men 
alle met digitale hulpmiddelen uitgevoerde bewerkingen 
op beelden rekenen. Hierbij wordt het begrip 'beeld' 
zeer ruim opgevat als twee-dimensionale data. Voorbeel 
den hiervan zijn fotografische opnamen van met het blote 
oog waarneembare verschijnselen, maar ook opnamen die 
verkregen zijn met behulp van röntgen- of gammastraling, 
radar of sonar, elektronenmicrescopie, radio-astronomie, 
etc. Met betrekking tot het doel van de verwerking onder- 
scheidt men in her algemeen (Rosenfeld 1976a, Gonzalez, 
1977, Pratt 1978, Castleman 1979) een viertal deelgebie- 
den binnen de digitale beeldverwerking: beeldeodering 
(voor een zo efficient mogelijke transmissie of opslag 
van beelden), bee!dreeonestructie (het terugwinnen van 
een beeld vanuit een vervormde en/of verruiste versie op 
grond van mathematische modellen), beeldverbetering (op 
empirie gebaseerde technieken om de voor een bepaald pro- 
bleem relevante beelddelen te benadrukken) en beeldana- 
lyse (analyse van een beeld in termen van de door het 
beeld gerepresenteerde verschijnselen). Een belangrijke 
stap in het beeldanalyse-proces wordt gevormd door de 


beeldsegmentatie. 


BEELDSEGMENTATIE 


Bij het coderen, reconstrueren of verbeteren van beelden 
wordt het resultaat van de bewerkingen opnieuw gevormd 
door een beeld. Dit is bij beeldanalyse niet het geval. 
De bewerkingen resulteren nu in uitspraken omtrent de fe- 
nomenen die door het beeld gerepresenteerd worden. Als 
we bijvoorbeeld röntgenopnamen van de kransslagaders van 
het hart analyseren, kunnen de resultaten van de digitale 
analyse bestaan uit uitspraken omtrent de positie en 
ernst van vernauwingen in deze kransslagaders (Gerbrands 
1980), Om dergelijke uitspraken te kunnen doen is het 
noodzakelijk een beschrijving van het beeld op te stel- 
len in termen van de verschillende gebieden of componen= 
ten waarin het beeld ontleed kan worden. In het beeld- 


analyse-proces kan dan ook veelal een eerste stap van 


Fig. 


beeldsegmentatie worden onderscheiden. Bij deze segmen- 
tatie wordt van ieder beeldpunt vastgesteld tot welk 
deelgebied het behoort. Indien deze gegevens in een ma- 
trix worden geplaatst, vormt deze een plattegrond van 
het beeld waarop de deelgebieden in kaart zijn gebracht. 
In de verdere analyse wordt van deze plattegrond ge- 
bruik gemaakt om de gebieden te classificeren of te be- 
schrijven. In de beeldsegmentatie=stap blijven de gebie- 
den onbenoemde en betekenisloze elementen, waaraan pas 
op hogere niveaus in het veelal hiërarchisch opgebouwde 
analyse-proces betekenis wordt toegekend. De plaats van 
de segmentatie-stap in het analyse-proces wordt geïllus- 


treerd in het sterk vereenvoudigde schema van Fig. Ì. 


beschrijving 


analyse object-relaties 


Pae, 


Plaats van Beeldsegmentatie binnen beeldanalyse, 
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Het is natuurlijk wel de bedoeling dat het segmentatie= 
resultaat interpretatie op hogere niveaus mogelijk maakt. 
Dit betekent dat de ontwerper doelgericht en probleemge- 
bonden te werk moet gaan bij het inrichten van de segmen=- 
tatieprocedure. 

Los van een specifieke probleemstelling wordt de 
doelstelling van beeldsegmentatie geformuleerd als het 
verdelen van een beeld in deelgebieden die "homogeen zijn 
in een bepaalde eigenschap". Elk beeldpunt a, uit de ver 
zameling beeldpunten A = { As As 5 an! krijgt één 
uniek etiket of label A, uit de verzameling labels A met 
Am { A, A e A} toegewezen dat aangeeft tot welk 
deelgebied het beeldpunt behoort. Er bestaat een grote 
verscheidenheid aan beeldsegmentatie=technieken om deze 
label-toekenning te verzorgen. De meeste hiervan kunnen 
worden gecombineerd met een relaxatie-procedure met als 
doel contextuele inconsistenties in de label-toekenning 


te reduceren. 
RELAKATIE-LABELPROCEDURE 


Relaxatie-labelprocedures zijn iteratieve procedures die 
er op gericht zijn contextuele informatie te gebruiken 
om lokale inconsistenties in de label-toekenning te redu- 
ceren. We gaan hierbij uit van een initiele label-toeken- 
ning, waarmee bedoeld wordt dat aan ieder beeldpunt a,EA 
een verzameling mogelijke labels A, c A wordt toegekend. 
Voorts dienen we te beschikken over kennis omtrent de 
verenigbaarheid of compatibiliteit van elk label van een 
beeldpunt a, met de labels van de buurpunten van a, 

Een label-toekenning heet ondubbelzinnig als aan elk 
beeldpunt a, precies êën label is toegekend, dat compati- 
bel is met het label van elk van de buren. Uitgaande van 
een dubbelzinnige initiële label-toekenning met meer dan 
Eén mogelijk label per beeldpunt worden nu de labels van 
naburige punten in rekening gebracht via de compatibili- 
teitsrelaties teneinde de dubbelzinnigheid te verminderen 
of geheel weg te nemen. De initiele label-toekenning 
wordt verkregen met behulp van operatoren die werken op 
het niveau van de beeldpunten van het aangeboden beeld en 
is dus data-afhankelijk. De compatibiliteits-relaties 
drukken de verenigbaarheid van labels uit. Ze worden in 
het algemeen gekozen voor de gehele klasse van te segmen- 
teren beelden en bieden de gelegenheid onze heuristische 
a priori-kennis omtrent het probleem in de procedure in 
te brengen. 

In de zogenaamde discrete relaxatie-methoden is een 
label voor een bepaald beeldpunt ôf mogelijk öf onmogelijk 
en zijn labels voor naburige beeldpunten ôf verenigbaar 
óf onverenigbaar. Het discrete relaxatie-algoritme verwij- 
dert dan uit de initiële label-verzameling van elk punt 
de labels die niet met ten minste êên label van elk van 
de buren verenigbaar zijn, om zo te komen tot een ondub- 
belzinnige label-toekenning (Rosenfeld 1976b). 


We beperken ons nu verder tot de zogenaamde probabi- 


listische relaxatie-methoden (Rosenfeld 1976b, Zucker 
1978). Hierbij wordt een label A, in zekere mate of met 
een zekere betrouwbaarheid aan een beeldpunt a, toegekend, 
hetgeen wordt uitgedrukt in weegfactoren P. (A). Deze 
weegfactoren voldoen aan 

PA)= 0 (1) 

EPA = | (2) 
Hierdoor kan men PA) opvatten als een schatting van de 
kans dat label A, het correcte label is voor beeldpunt a, 
In een probabilistische relaxatie-procedure worden nu deze 
schattingen iteratief bijgesteld op grond van de schattin- 
gen van de kansen van de labels van de naburige punten en 
de bijbehorende compatibiliteits-coëfficiënten. De compati- 
biliteits-coëfficient van label A, voor beeldpunt a, en 
label en voor beeldpunt ä, wordt aangeduid met Cz) 
Deze compatibiliteits-coefficienten mogen waarden aannemen 
in het interval [-1,1l. Hierbij betekent CA) =| 
dat label A, voor a, zeer verenigbaar is met label A, voor 
beeldpunt a, een waarde -l betekent dat label A, voor a, 
onverenigbaar is met A, voor a, terwijl Cs Are) = 0 
aangeeft dat label A, voor a, niet beïnvloed wordt door \, 
voor dr Voor overige situaties kunnen tussenliggende 
waarden worden gekozen. 

Bij het iteratief bijstellen van de schattingen van 
de kansen wordt gebruik gemaakt van een bijstel-regel van 
de vorm 

PL OP = B PAD s Ars ) (3) 
waarin t het nummer van de iteratie aangeeft en Ai de 
invloed van de omgeving. In iedere iteratie worden de 
schattingen voor alle labels A, voor alle beeldpunten a, 
bijgesteld met een in principe parallel algoritme. Hier 
betekent parallel dat binnen een iteratie de resultaten in 
een bepaald punt niet beïnvloed worden door de binnen die 
iteratie elders verkregen resultaten, zodat alle bewerkin- 
gen simultaan zouden kunnen worden uitgevoerd. Een veel 
gekozen mogelijkheid voor de functie ® { s … } is de 
volgende niet-lineaire bijstel-regel: 

PODIA) 


En) (4) 


u : PAD [A+ A | 
waarin de omgevingsinvloed wordt gedefinieerd als 
bn Fay Eje) PAP 6) 
De sommatie over j in (5) wordt uitgevoerd over de beeld 
punten 4, in de te beschouwen omgeving van het beeldpunt 
a,. Met de coëfficiënten d,. kan de totale invloed die a. 
heeft op a, worden gewagen, Kiest men deze weegcoëfficien- 


ten zodanig dat geldt 


NE (6) 
J. 4 
dan volgt daaruit dat 


waardoor te allen tijde voldaan is aan de eis (1) dat de 
schattingen van de kansen niet-negatief zijn. Daarnaast 
draagt de noemer in (4) zorg voor de eis (2) dat de som 
van de schattingen voor de kansen van de labels in ieder 


beeldpunt steeds gelijk blijft aan |. Eigenschappen als 


convergentie en stabiliteit van probabilistische relaxa- 
M(i) 


tie met behulp vaa bijstel-regel (4) worden beschreven P,(R) = MAX Mu) (11) 
in (Zucker 1978). In de volgende paragrafen behandelen KL 
we een tweetal voorbeelden. waarbij het maximum in de noemer wordt bepaald over het 


gehele beeld, De schatting van de kans dat a, gêën rand- 


VOORBEELDEN punt is wordt dan 
Rand-dectectie P‚(R) = LP, (R) (12) 
Indien de in het beeld te onderscheiden gebieden van el- De compatibiliteits-coefficienten dienen nu zodanig geko- 
kaar verschillen in (gemiddelde) helderheid, kan gebruik zen te worden dat PR) wordt verhoogd als een randelement 
gemaakt worden van segmentatie-methoden die gebaseerd in a, ook wat richting betreft ondersteuning vindt in de 
zijn op het detecteren van de randen van de gebieden. De- omgeving en juist verlaagd wordt als dat niet het geval is. 
ze randen komen overeen met veranderingen in de helder- De richting van het randelement door een punt wordt lood=- 
heid, die in principe met differentiërende operatoren recht genomen op de richting van de gradiënt in dat punt. 
gedetecteerd kunnen worden. Het verdient aanbeveling deze De richting van de rand in het te beschouwen punt a, wordt 
operatoren zo in te richten dat de ruisinvloedenenigs- aangeduid met @, de richting in een buurpunt ej met B, de 
zins beperkt warden. In dit voorbeeld maken we dan ook richting van de verbindingslijn van a, naar ä, met Y en 
gebruik van de volgende operatoren B, en Ö de afstand tussen a, en a met d. Door Schachter c.s, 
(Schachter 1977) zijn nu de volgende compatibiliteits- 
+1 +l +1) +1 0 —lI coëfficiënten voorgesteld: 
6‚=zjo o of s-jn 0 (8) N 
el ne +1 0 =I G,;(R‚R): cos (ar) .cos (Br) .2 (13) 

Het oorspronkelijke beeld wordt twee-dimensionaal gecon- CBD: min[0, cos (2u 2) 2 (14) 
volueerd met beide operatoren, waarbij het resultaat 
6, (i) respectievelijk ö(i) wordt toegekend aan het beeld- C,; RD): (l-cos (26-27) „2 (À*DD (15) 
punt a. waarvan de positie overeenkomt met het centrum 
van de operator. De operator zn heeft een differentië- cj): ze (16) 
rend karakter in vertikale richting en is dan ook het 
meest gevoelig voor horizontale randen, De operator Ö, Het effect van deze keuze van de compatibiliteits-coeffi- 
is juist het meest gevoelig voor vertikale randen. Beide ciënten wordt geïllustreerd met een aantal voorbeelden in 
operatoren samen kunnen worden gecombineerd in de rich- Fig. 2. In deze figuur komt het gesloten pijltje overeen 
tingsongevoelige gradient waarvan de modulus M(i) en de met het randelement door het te beschouwen punt a;, 
richting 8(i) worden gegeven door: het open pijltje met de rand door a. terwijl het label R 

M(i) = A6, Li} + 16,1} (9) 

y x 
8 (1) 
Bli) = arctan 5 TD (10) 
x 

Als een beeldpunt a, op de rand tussen twee gebieden ligt, 
is de waarde van M(i) relatief hoog. Als randpunten wor- 
den dan ook gewoonlijk die punten a, gedetecteerd waar- 
voor geldt dat M(i) een bepaalde drempelwaarde over- 
schrijdt. In termen van het toekennen van labels aan de 
beeldpunten komt dit overeen met het gebruik van twee 
labels, hier aangeduid met R(rand) en R (geen rand), die 
worden toegewezen op grond van de modulus van de gra- 
dient. Teneinde nu in deze procedure van rand-detectie 
een probabilistische relaxatie-fase op te nemen, kunnen 
de initiele schattingen van de kansen op de labels R en 
R worden gebaseerd op de gemeten waarden van de modulus Fig. 2: Voorbeelden van de compatibiliteits-coëfficiënten 


van de gradiënt volgens 


wordt aangeduid met een st LD. Een positieve waarde van 


de compatibiliteits-coefficient betekent versterking, 


een nesatieve verzwakking van het betreffende label. 


Fig. 4: Modulus van Fig. 5: Relaxatie-resul- 


gradient taat (5 iteraties) 


Ter 1llustratie passen Wet de hier reschets te methode toe 
op het testheeld van Fig. 3. De modulus van de gradient 
Is gegeven in Fig. & en het relaxatie=resultaat na vijf 
iteraties in Fig. 5. In beide figuren is de helderheid 

in een heeldpunt evenredig met de schatting van de kans 


dat het beeldpunt een randpunt is. Fig. 6 representeert 


omgeving in de relaxatie=procedure wordt gevormd door de 
acht buren van elk beeldpunt. Opgemerkt moet nog worden 
dat ook de richting f(i) van de gradiënt iteratief werd 
bijgesteld en daarmee de hoeken @ en B van de randelemen=- 


ten. Voor details wordt verwezen naar (Schachter 1977). 


Cluster=segmentatie 
[In dit tweede voorbeeld van het gebruik van relaxatie- 
procedures in de beeldsegmentatie gaat het om een segmen- 
tatie-methode waarbij alle beeldpunten worden geclassifi- 
ceerd in een aantal klassen met behulp van technieken 
uit de statistische patroonherkenning (Duda 1973). Hier- 
bij wordt van elk beeldpunt een aantal kenmerken bepaald, 
die worden geplaatst in een kenmerkvector. Door middel 
van deze kenmerkvectaor wordt ieder beel dpunt afgebeeld 

in de kenmerkruimte. Afhankelijk van zijn positie in de 
kenmerkruimte wordt vervolgens ieder beeldpunt geclassi- 
ficeerd tot een van de te onderscheiden klassen. 

Indten men geen voorkennis heeft over de omliggende 
statistische verdelingen in de kenmerkruimte, noch de 
beschikking heeft over een voldoend groot aantal a priori 
gelabelde leersamples, kan men gebruik maken van zoge- 
naamde cluster-methoden bij het ontwerpen van het classi 
ficatie=-mechanisme. Hierbij wordt uitgegaan van de ken= 
merkvectoren van een verzameling ongelabelde beeldpunten, 
die worden afgebeeld in de kenmerkruimte. Met behulp van 
een cluster-algorithme brengt men nu een partitie teweeg 
van de verzameling in een aantal homogene clusters. De 
term homogeen betekent hier dat beeldpunten in hetzelfde 
cluster in hun kenmerkbeschrijving onderling meer gelij 
kenis vertonen dan ten opzichte van punten iù andere 
clusters, Van de gevonden clusters worden de centra be- 
paald, die vervolgens bijvoorbeeld in een minimum-af- 
standsclassifticator gebruikt kunnen worden. 

In dit voorbeeld wordt als eerste kenmerk de helder- 
heid of grijswaarde van de beeldpunten genomen. Als twee= 
de kenmerk kiezen we een zogenaamde textuurmaat. Textuur 
kan worden omschreven als de lokale variaties in de grijs- 
niveaus, die karakteristiek zijn voor de aard van een 
waarsenomen oppervlak, Textuurmaten zijn er op gericht 
deze karakteristieke lokale variaties te kwantificeren. 
Er bestaat een groot aantal van deze textdurmaten 
(Haralick 1979), In dit voorbeeld gebruiken we de door 
Haralick voorgestelde maat GLD-ASM om in een 5x5 venster 
rond ieder beeldpunt de textuur te kwantificeren, waarna 
de gevonden textuurwaarde steeds als kenmerk wordt toe- 
gevoegd aan het beeldpunt dat overeenkomt met het centrum 
van het venster. Op deze manier wordt aan alle 65536 
beeldpunten van het testbeeld van Fig, 3 een twee=dimen= 
sionale kenmerkvector toegevoe ed, 

Met een aselecte steeknruef wan 6500 heel dpunten 15 
vervolgens de cluster=analyse uitgevoerd. Hierbij is ge- 
bruik gemaakt van het zogenaamde "k=-means cluster algo- 


rithm”’ (Ball 1967). Het principe van dit algorithme is 


als volgt. Kies k (bijvoorbeeld 5) initiele clustercen- 
tra in de kenmerkruimte. Vervolgens worden de data-punten 
toegewezen aan een clustercentrum op grond van een mini= 
mum-afstandscriterium. Van de op deze wijze gevonden 
clusters worden weer de centra bepaald. De procedure 
wordt herhaald tot de clustercentra niet meer verschui= 
ven. Toepassing van dit cluster-algorithme leidde in ons 
geval tot drie clusters. 

Classificatie van alle 65536 beeldpunten op grond van 
een minimum-afstandscriterium ten opzichte van de drie 
clustercentra resulteert in het in Fig. 8 weergegeven 
gesegmenteerde beeld. De segmentatie is nogal gefragmen- 
teerd. Het cluster=algorithme zoekt homogene clusters in 
de kenmerkruimte, maar dit behoeft nog niet te leiden 
tot een segmentatie in homogene gebieden in het beeld 
domein. De lokale consistentie in de label-toewijzing 
kan weer worden verbeterd door een probabilistische 
relaxatie-procedure. 

De initiële schattingen van de kansen op de verschil- 
lende labels worden afgeleid uit de afstanden van een 
datapunt in de kenmerkruimte tot de verschillende clus- 


tercentra volgens 


(17) 


waarbij c het aantal clusters is en d, (A) de afstand van 
het met beeldpunt a, overeenkomende data-punt in de ken- 
merkruimte tot het centrum van het met label A, aange= 
geven cluster. 

De compatibilsteits-coefficienten worden hier geko- 
zen vanuit de gedachte dat we streven naar relatief 
grote gebieden met eenzelfde label. Vanuit die achter- 


grond maken we in dit voorbeeld gebruik van 


|l als \, mt A. 
Cz Oro) = (18) 
=| als A, 4 A, 


omdat het meer voor de hand ligt dat buurnunten in het 
beeld hetzelfde label hebben dan een verschillend label. 
Het is duidelijk dat deze keuze van de compatibiliteits- 


coëfficiënten weinig subtiel is, waardoor o.m. langge- 


Fig. 8: Clusterresultaat Fig. 9: Clusteren gevolgd 
door relaxatie 


overlevingskans hebben. 

Het resultaat van de probabilistische relaxatie-label- 
procedure met de in (17) gedefinieerde initiele kansen 
en de in (18) gegeven compatibiliteitscoefficienten is 


weergegeven in Fig. 9 
SLOTOPMERKINGEN 


rekte smalle gebieden (de statiefpoten in Fig. 3) weinig 

Uit de hier gegeven voorbeelden moge blijken dat probabi- 

listische relaxatie-methoden zeker van nut kunnen zijn bij 

het segmenteren van beelden. In feite wordt hiermede de 

segmentatie=procedure een hiërarchische procedure waarbij 

op het niveau van de oorspronkelijke beeldpunten de ini- 

tiele schattingen van de labelkansen worden opgesteld, | 

terwijl op een hoger niveau deze schattingen worden bij- | 

gesteld op basis van de lokale context. Opgemerkt dient 

nog te worden dat relaxatie-procedures rekenintensief 

zijn. De rekentijden liggen in de orde van | minuut per | 

iteratie voor een 256x256 beeld met de hier gehanteerde | 

experimentele FORTRAN-programma's op een IBM-370/158 | 

rekensysteem. | 
De auteur hecht eraan, de samenwerking met 


P.C.W. Sommen, H.A.J. Merkus, R.M. Bolle, E. Backer en 


IJ. Boxma op deze plaats te noemen. 
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OPTIMAAL STATISTISCH FILTEREN VAN BEELDEN 


Ir. J. Biemond 


Technische Hogeschool Delft 


Afdeling der Elektrotechniek, Vakgroep Informatietheorie 


Optimal 


statistical filtering of images, An optimal line-by-line recursive Kalmanfilter is derived for 


restoring images degraded by spatially invariant degradation phenomena in the presence of additive white 


noise, Results of several filter experiments on noisy motion blurred images and on noisy defocussed 


images are presented. 


INLEIDING 

Afbeeldingen of beelden van een scène (in de optica: ob- 
ject) zijn in meer of mindere mate behept met vervorming 
en ruis, Dit komt omdat geen enkel beeldvormend systeem 
(elektronenmicroscoop, camera, telescoop) een perfect, 
onvervormd beeld afgeeft. Ruis kan geïntroduceerd worden 
tijdens de opname, in de transmissieweg of gedurende het 
registratieproces (korrelruis), terwijl bij vervorming 
gedacht kan worden aan sferische aberratie van lenzen, 
bewegingsonscherpte, onjuiste focussering, niet-lineaire 
eigenschappen van fotografische materialen, geometrische 
vervorming of — bij telescoopfoto's — vervorming ten ge- 
volge van atmosferische turbulenties. 

Is een afbeelding van groot belang om beslissingen 
te nemen, of erg zeldzaam, (satellietbeelden, elektronen- 
microscoopbeelden, historische opnamen), dan kan gevraagd 
worden naar technieken die de vervorming en de ruis mini 
maliseren en een gereconstrueerd beeld opleveren, Men de- 
finieert beeldreconstructie dan als een vorm van (digi 
tale) beeldverwerking, waarbij men het oorspronkelijke 
beeld "zo goed mogelijk", d.w.z, gebruik makend van een 
objectief evaluatiecriterium, tracht terug te winnen uit 
een vervormde en door ruis gestoorde versie, Dit is een 
bewerking achteraf, men gaat er vanuit dat het beeld- 
vormend systeem een gegeven 18, waar men geen invloed 
meer op heeft, 

In figuur | is het hier te hanteren model van het 


beeldvormende systeem geschetst: 


vii,j) 


Fig. |. Beeld(ver)vormend systeem. 


Het origineel x(i,j) dat in de optica aangeduid wordt als 


het object, doorloopt een twee-dimensionaal (2-D) beeld- 


vormend systeem. Het uitgangsbeeld y(i,j) is het beeld dat 
wij ter beschikking hebben. De bedoeling is nu, uitgaande 
van yli,j), het origineel x(i,j) te reconstrueren. We be- 
perken ons tot monochrome, discrete beelden, die gerepre- 
senteerd worden door het helderheidsniveau x(i,j) (resp. 
yv(i,j)) in ieder punt (1,j) met i de rij-index en j de 
kolomindex in het gebied 1 Si SM, 1Sj SN, 

Het is belangrijk om op te merken dat beeldrecon- 
structie a priori kennis vereist omtrent het originele 
beeld, de aard van de vervorming en de ruis. Zonder enige 
a priori kennis is het onmogelijk om zelfs maar te ont- 
dekken dat er iets fout is met een gegeven beeld, Het is 
immers de a priori kennis die ons leert dat de korrel 
structuur van een foto een artefact is en dus eigenlijk 
niet op de foto thuishoort. 

Omtrent het originele beeld zouden we graag structu- 
rele informatie bezitten, De werkelijkheid, zoals die be- 
sloten ligt in een beeld is meestal zo complex, dat a 
priori deterministische of linguistische modellen buiten 
bereik blijven. Daarom hanteert men vaak eenvoudige sta- 
tistische beschrijvingswijzen: het onvervormde beeld wordt 
gemodelleerd als het 2-D scalaire uitgangssignaal van een 


lineair systeem L,, gestuurd door Gaussisch witte ruis, 


; 

De verona: die door het beeldvormend systeem 
geïntroduceerd wordt, is in het algemeen zeer complex. 
Echter in veel gevallen van praktische betekenis, zoals 
lineaire bewegingsonscherpte, defocussering en atmosferi- 
sche turbulenties kan een lineaire systeembeschrijving ge- 
hanteerd worden, Een vervormd beeld kan dan worden ver- 
kregen als uitgangssignaal van een lineair systeem Ls 
met als ingangssignaal het onvervormde beeld, 

De observatieruis is meestal stochastisch van aard, 
en wordt beschreven met behulp van een kansdichtheid. We 
veronderstellen de ruis additief Gaussisch wit, onafhan=- 
kelijk van het beeld, en toegevoegd aan het reeds ver- 
vormde beeld. 

Zodoende kan een vervormd en door ruis gestoord beeld 
yÜi,j)gemodelleerd worden als het uitgangssignaal van een 


lineair systeem L (L, en L, in cascade), gestuurd door 


Ì 
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witte ruis, waaraan vervolgens additieve observatieruis 
wordt toegevoegd (zie fig, 2). 


Observ. ruis 


L 
mm en mn ET 
mn 


vii,j) [Lineair Lineair beeld- 
systeem L, vormend syst. L 


k Onvervormd J 
| beeld | 


Fig. 2: Vervormd en door ruis gestoord beeld, gemodel- 


x(i,j) 


leerd als uitgangssignaal van een lineair sys- 


teem gestuurd door witte ruis, 


Vanwege de statistische beschrijving van zowel 
beeld als ruis wordt het reconstructieprobleem opgevat 
als een optimaal filterprobleem, Door de aanwezigheid van 
observatieruis zijn inverse filtermethoden vrijwel on- 
bruikbaar, Het gereconstrueerde beeld wordt verkregen 
aan de uitgang van een lineair filter, waaraan de ver- 
vormde en verruiste observaties worden toegevoerd, Het 
filter wordt hierbij zodanig gekozen, dat het gemiddelde 
kwadratische verschil tussen het oorspronkelijke beeld 
x(1,j) en het gereconstrueerde beeld %(i,j): 


e“ m Bl{xCE, DR, 1932 1) 


geminimaliseerd wordt (zie fig, 3). 


Optimaal 
lineair filter 


Ri, j) 


Waargenomen beeld gereconstrueerd 


beeld 


Fig. 3: Beeldreconstructie als optimaal lineair filter- 


probleem, 


Een tweetal optimaal statistische filters komt in aan- 
merking: het Wienerfilter en het Kalmanfilter. Beide fil- 
ters gebruiken dezelfde a priori kennis, n.l. de auto- 
correlatiefuncties van beeld en ruis, en de overdrachts- 
functie van het beeld(ver)vormende systeem, Bovendien 
minimaliseren ze beide de gemiddelde kwadratische fout. 
Zodoende kan dezelfde reconstructiekwaliteit verwacht 
worden, Toch is gekozen voor het Kalmanfilter om de 
volgende redenen: 
(1) het Kalmanfilter is recursief implementeerbaar, 
(2) de toepassing van het Kalmanfilter is in principe 
niet beperkt tot stationaire processen. 
Vaor de afleiding van het Kalmanfilter is een line- 
air dynamisch model nodig, dat enerzijds een recursieve 


beschrijving geeft van de statistische eigenschappen van 
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het oorspronkelijke beeld (het beeldmodel L) en dat an- 
derzijds de vervorming en de ruis modelleert (het beeld 
vormend systeem LL). We zullen achtereenvolgens op beide 


modellen L, en L, wat dieper ingaan. 


RECURSIEVE BESCHRIJVING VAN HET BEELDMODEL Li 
Ee ande 
Door een beeld op te vatten als een realisatie van een 


homogeen discreet stochastisch veld 
TEE II Telde nere HEL Does AT 


is het mogelijk de eerste- en tweede-orde eigenschappen 
van dit veld te gebruiken als a priori statistische kennis 


voor het ontwerpen van recursieve beeldmodellen: 


Elx(i,j)} = u 
(2) 
R(k‚1) = Elx(i,j)x(itk, j+1)] 


Twee veelvuldig in de literatuur gebruikte autocorrelatie- 
functies om een klasse van beelden te beschrijven zijn 


(Jain, 1977): 
(1) Scheidbaar negatief-exponentiële autocorrelatiefunctie: 


Rík,1l) = a D, 


2, ll Tl, on 


2 Il Ld 
met JO de variantie van het beeld en p, en Dn de ver- 


tikale en horizontale correlatiecoefficient (zie fig. 4). 


(2) Isotrope correlatiefunctie 


Alet 


Rík,1) = p (4) 


met p decorrelatiecoëfficient en s een vormparameter (zie 
fig. 5). 


Fig. 4: Scheidbare negatief- exponentiële autocorrelatie- 


functie, 


Fig. 5: Isotrope autocorrelatiefunctie s=l, 


De vraag dient zich nu aan, hoe op basis van de statis- 
tische gegevens in (2) een recursief beeldmodel opgesteld 
kan worden in de vorm van een stochastische differentie= 
vergelijking. In analogie met het 1-D geval zou gedacht 
kunnen worden aan het schema in fig. 6. Dit 

schema is echter niet in zijn geheel bruikbaar door het 
ontbreken van een algemeen 2-D spectraal factorisatie- 
theorema, Concluderend kan gezegd worden dat het uit- 
gaande van een geschikte eerste- en tweede-orde beschrij= 
ving in het algeneen niet mogelijk is één-éénduidig 

een stochastische differentievergelijking op te stellen. 
Het probleem wordt in de literatuur daarom meestal om- 
gekeerd, Uitgegaan wordt van een algemene stochastische 
differentievergelijking als beschrijving van een 2-D 
stochastisch veld: 


x(i,j) = EF alp,g)x(itp,jta)tvli,j) (p,ES (5) 


PsJ 

d.w.z. de helderheid van een beeldpunt x(i,j) wordt ge- 
schreven als een eindige som van omliggende beeldhelder- 
heden x(itp,j+q) en een stochastische term v(i,j). Deze 
stochastische term kan gezien worden als ingangssignaal 
van het proces of als de Fout die gemaakt wordt bij het 
genereren van x(i,j) met behulp van x(itp,j+q) met 
(p‚a) € 8. S is zen coöordinatenverzameling met uitslui 
ting van (p,q) = (0,0), Met de keuze van de verzameling 
S is de structuur van het beeldmodel vastgelegd. Aange- 
nomen wordt dat «(i,j) en vli,j) gemiddelde waarden nul 
bezitten, Door een geschikte keuze van de coefficienten 
alp,q) kan men een model zo goed mogelijk aanpassen aan 


de gegeven eerste= en tweede-orde statistiek, bijv. door 


de coefficienten zodanig te kiezen dat de variantie van 
e(i,j) geminimaliseerd wordt, Welke representant van het 
algemene beeldmodel uiteindelijk als beeldmodel t.b.v, 
het Kalmanfilter zal dienen hangt af van de rekentechni- 
sche implementatie en de kwaliteit van het gekozen model. 
Om de kwaliteit te meten zijn een aantal criteria ontwik= 
keld, waaronder de Cramêr-Rao grens (Biemond, c.s., 1979). 
Een aantal in de literatuur gebruikte beeldmodellen 


wordt hier gegeven: 


(1) Lijnseanningsmodel (Nahi, Assefi, 1972) 
jr 3 


z(i,j) = axli,j-l)sbvli,j=1) (6) 


(2) DPCM-model (Habibi, 1972) 


Je 
i=l Db @ x(i,j) = ajx(i-l,j)ta,x(i,j-I)+ 
1 GO & e= ajagx(iel, j-l)+bv(iel, j-1) (7) 


(3) Semi-ecausaal model (Jain, 1977) 


dek 4 


il. @ @ 
1. 0 
El. Ei OO 


x(i,j) = a,{x(i-l, j)exli+l, jh 
+ a,lxli-l,j-Itx(i+l, je 
+ agx(i,j-l)tbv(i, j) (8) 


(4) Niet-eausaal nearest neighbour model (Jain, 1974) 


sel j 11 


iel. « @ « _« xli,j) = a, {xli-l,j)txliel, jj 
Ez 8 & & Ò + a,lxli,j-Dtrli,j+D bv{i,j) 
BEL. Oan OR a (9) 


Uit de naamgeving en de structuur van de modellen 
valt af te leiden dat we ons bedienen van een of andere 


definitie van causaliteit in 2-D, We zullen daarom kort 


KL} + RI) > Sluv) > Hjusiv) => hlk‚D) > x(i,j) = EF alp‚g)x(itp,j+g)tv(i,j) 
stochastisch autocorre- | vermogens- overdrachts | impuls PA stochastische 
veld latie diechtheids- | functie responsie differentie vergelijking 
functie spectrum 
Fourrier spectrale inverse 
transfor- factorisatie Fourrier 
matie transfor- 
matie 


Fig. 6: Spectraal factorisatie schema in 2-D, 
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ingaan op de betekenis van causaliteit in 2-D en een 
tweetal definities geven. In l-D betekent causaliteit 
dat een signaalwaarde slechts een functie is van voor- 
gaande signaalwaarden en een ingangsterm. In 2-D daar- 
entegen is er geen natuurlijke ordening van punten, zo- 
dat een 2-D vlak oneindig veel generatierichtingen door 
een punt x kan bevatten. Bovendien is niet duidelijk wat 
heden, verleden en toekomst binnen een beeld is, omdat 
geen van beide onafhankelijke variabelen de tijd is. Om- 
dat recursieve filteralgorithmen, gebaseerd op dit type 
modellen, in het tijddomein werken, moet de tijd op 
kunstmatige wijze geïntroduceerd worden. Twee definities 
van causaliteit worden gehanteerd: 
— Non-symmetric half-plane (NSHP) causaliteit (Woods- 
Radewan, 1977). 
Dit type causaliteit is gebaseerd op het onder 1) ge- 
noemde lijnscanningsmodel, Alle beeldlijnen worden 
achter elkaar geplaatst, zodat het beeld op de volgen- 
de wijze in een heden, verleden en toekomst wordt ver- 
deeld (zie fig. 7). 


Verleden 


Fig. 7: NSHP-causaliteit 


—- Quarter-plane (QP) causaliteit (Habibi, 1972), 
Dit type causaliteit is gebaseerd op het onder 2) ge- 


noemde DPCM-model (zie fig. 8). 


Toekomst 


Fig. 8: QP-causaliteit 


HET BEELD(VER)VORMENDE SYSTEEM La 


De vervorming t.g.v. het beeldvormend systeem be- 
schrijven we met behulp van een lineair systeem. Aan de 
ingang van dit systeem wordt het originele onvervormde 
beeld aangeboden, terwijl aan de uitgang het vervormde 
beeld aanwezig is. De beeldruis die op verschillende 
plaatsen in het beeldvormende systeem ontstaan kan zijn, 
wordt gemodelleerd als een additief signaal-onafhankelijk 


proces, dat nä de vervorming optreedt. De lineaire beeld- 
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vorming wordt dan gemodelleerd m.b.v. de volgende 2-D 


superpositiesom: 
sjd= E ax(k, Dali: k,‚1)tw(i, j) (10) 
k 


met y(i,j) het waargenomen beeld op het punt (i,j), x(k,1) 
het originele beeld op het punt (k,1), wli,j) de signaal- 
onafhankelijke ruis en gli,j; k,1) de impulsresponsie of 
puntspreidingsfunctie, die het uitgangsbeeld beschrijft 
bij een puntbron aan de ingang. Is het uitgangsbeeld, af- 
gezien van een translatie, onafhankelijk van de positie 
van de puntbron aan de ingang, dan heet het systeem ver- 
schuivingsinvariant en reduceert (10) tot een convolutie- 


SO ; 


yi, = & elk, Dgldik, j-1)tw(i „j) 
CHI) 
= k Ex(i-k,j-1)g(k, 1)+w(i, j) 
k 1 
Voorbeelden van verschuivingsinvariante vervormingen zijn 
de eenparige beweging langs een as en de defocussering. 
Een puntbron aan de ingang geeft dan een streepje resp. 
een vlek in de opname y(i,j), onafhankelijk van de posi- 
tie van de puntbron aan de ingang. Heeft men a priori 
kennis omtrent de oorzaak van de vervorming, dan is het 
mogelijk de puntspreidingsfunctie door berekening te be- 
palen, waarna sommige parameters, zoals de afstand waar- 
over de beweging heeft plaats gevonden òf de diameter van 
de defocussering uit het beeld y(i,j) worden geschat. 


Ter illustratie dient het volgende voorbeeld. 


Voorbeeld 


Van de oorspronkelijke helderheidsverdeling (scène) 
x(i,j) wordt een fotografische opname gemaakt, waarbij de 
camera met constante snelheid V wordt bewogen in de j- 
richting. Deverkregen opname vertoont bewegingsonscherpte. 
Als de belichtingstijd bekend is en wordt gegeven door 
5 


het interval 55), dan geldt in het continue geval voor 


de puntspreidingsfunctie (Gerbrands, Backer, 1980) 


Blikje 7 -WTS jl SWT; ik = 0 
0 


(12) 
elders 

In figuur 9 is deze puntspreidingsfunctie getekend. Is 

de belichtingstijd onbekend, dan kan deze achterhaald 

worden door het gedrag van de vervorming te analyseren 

in het Fourierdomein, Fouriertransformatie van (12) 

geeft de overdrachtsfunctie G(u,v): 


sin 2VVT 


Glu,v) = (13) 


buVT 
Dit is een sincfunctie in de v=-richting onafhankelijk van 
de u-richting. De modulus van de overdrachtsfunctie is 


geschetst in fig. 10, 


Fig. 9: Puntspreidingsfunctie van de lineaire bewegings- 


onscherpte langs de j-as. 


Fig. 10: Overdrachtsfunctie lineaire bewegingsonscherpte 


Het nulpuntenpatroon van deze overdrachtsfunctie bestaat 
uit evenwijdige rechte lijnen, welke loodrecht op de be- 
wegingsrichting staan. In het vermogensdichtheidsspectrum 
van het door lineaire beweginsonscherpte vervormde beeld 
manifesteren deze nulpunten zich als herkenbare minima. 
De afstand tussen deze vertikalen kan opgemeten worden 
en 1s direct gerelateerd aan de afstand waarover de ver- 
vorming zich heeft uitgestrekt. 

Opgemerkt dient te worden dat bij het verwerken van 
digitale beelden op een computer, gebruik gemaakt wordt 
van het discrete equivalent van de hier behandelde can- 


tinue puntspreidingsfunctie, 


DE KALMANF ILTER-/ERGELIJKINGEN 

De beide subsystemen L, en L, kunnen nu samengevoegd 
worden tot een dynamisch model, bestaande uit een beeld 
of toestandsvergelijking en een observatievergelijking. 
Op dit dynamisch model worden dan de Kalmanfilterverge- 
lijkingen geënt. 

Bij wijze van voorbeeld ontwerpen we een Kalmanfil- 
ter om lineaire bewegingsonscherpte en additieve signaal- 
onafhankelijke ruis uit een beeld te filteren. We kiezen 
het DPCM model om het originele onvervormde beeld te be- 


schrijven: 


KD > pri dep dl, -PpOyr j-De 


+ vis=l, jl) (14) 


De coefficienten po, en Pp, Zijn zodanig gekozen, dat dit 


h 


model een exacte beschrijving geeft van een klasse van 
beelden met een negatief-exponentiële scheidbare autocor= 
relatiefunctie. Modelvorming van de ruis en de lineaire 
beweging over de (discrete) afstand Io leidt tot de vol- 
gende observatievergelijking: 
Lo 
yi,j) = E_ x(i,j-Dgli)tw(i,j) (15) 
1=Û 

Pogingen om op dit 2-D dynamische model (14) en (15) een 
optimaal 2-D Kalmanfilter te baseren leiden schipbreuk. 
De resulterende schatters behouden niet de optimaliteit 
in de kleinste kwadratische zin van het l-D Kalmanfilter 
(Biemond, Gerbrands, 1980). Optimaliteit is wel gegaran- 
deerd als dit 2-D model wordt geschreven in matrix=vector= 
notatie (in êên onafhankelijke variabele) als een eerste= 
orde lineaire vector-differentievergelijking. Op basis 
van dit model kan een optimaal |-D vector-Kalmanfilter 
worden afgeleid dat lijn-voor-lijn filtert. 

Het opstellen van een matrix=vector-vergelijking 
voor het DPCM-model in (14) gaat als volgt. 


Eerst wordt (14) herschreven tot 
x(i, jp, x(i, jl) = PL riet, pxl, j- DH 
+ viel, j=1) (16) 


Stel vervolgens dat X(i) en V(i) kolomvectoren zijn met 
4 Lj Ll Ë tE - Pe 

KC) 5 Ebi, ldenli,2d, east, ND] @rl LEI Ee [v(i,l), 

wide es dan kan (16) in matrix-vector nota- 


tie geschreven worden als: 


AX) = p AXCi-ID+VCi-1) C17) 


met Â = 


De A-matrix is een beneden-driehoeks Toeplitz-matrix. De 
rijen zijn lineair onafhankelijk en dus bestaat de inver- 


se van A. Na vermenigvuldiging met a gaat (17) over in: 
WE 0 X(i-1+AT V(i-1) (18) 


Vergelijking (18), te zamen met geschikte beginvoorwaar- 
den, kan nu dienen als beeld- of toestandsvergelijking in 
het gebied{(i,j)met i=l,...M: j=l,....‚N}. 

De observatievergelijking (15) kan nu als volgt in 


vector-matrix notatie geschreven worden: 


Y(ij = GK(ij+Wij (19) 


ECI) 
et Gem p | ae 
me Saa Sa 
BUD Ez 
0 RCL) etl) 20) MxN 


Combinatie van (18) en (19) leidt tot het volgende dyna- 


mische model: 


Kil) 


p XLi)+AT Ivoi) 
(20) 


Y(i) GX ijeWji ) 


Hiervoor kunnen de volgende Kalmanfiltervergelijkingen 


worden afgeleid (Panda, Kak, 1977). 


X(i+1) = oC) 1) [Yr 1)-Go, KC) (21) 


predictie innovatie 


Uit deze vergelijking blijkt dat een schatting K(i+1) van 


de Ci+1)® beeldlijn wordt verkregen op basis van een 
schatting X(i) van de i° beeldlijn (predictie), waarna 
deze schatting gecorrigeerd wordt m.b.v. de nieuwe in- 
formatie in de observaties Y(i+l) op de Ci+1)° beeldlijn 
(de innovatieterm), gewogen met de versterkingsmatrix 
F{i+l), Voor deze versterkingsmatrix geldt de volgende 
relatie: 

P(iei) = QDEFLGAIIG EEJ! 
QCi+1) =p, AIF (i+1)G1Q(i)+AT 


AE en „It 
met Q(i)= pP(I)+A EA s 


1E (22) 


Ek 
Kd 


Eyehy en P(i) zijn de covariantiematrices van resp. 
Vi), Wi) en de foutvector [KI RC). De filterverge- 
lijkingen zijn bepaald onder de aanname dat het signaal 


een gemiddelde waarde nul bezit. 


EXPERIMENTELE RESULTATEN 


De prestaties van het bovenbeschreven vector-Kalmanfilter 


zijn experimenteel onderzocht voor een aantal S-bits 
testbeelden. Eerst is onderzocht in hoeverre het filter 
de ruis minimaliseert, Figuur Il is opgebouwd uit 4 
beelden ter grootte van báxb4 beeldpunten. Linksboven 
is afgebeeld het origineel (meisjesgezicht), rechts- 


boven hetzelfde testbeeld maar nu verontreinigd met door 


de computer gegenereerde additieve signaal-onafhankelijke 


Witte Gaussische ruis, met gemiddelde waarde nul en een 
signaal-ruisverhouding (SNR) van 2. De signaal-ruisver- 


houding is hierbij gedefinieerd als 


variantie signaal G: 


ONK ® Varlantie ruis Ee 
Figuur |l linksonder toont vervolgens het gefilterde 
beeld, gebaseerd op het DPCM-model met be 0,8, 
terwijl rechtsonder het gefilterde beeld is afgebeeld 
voor p, =p,” 0,9. Naast een goede ruisuitsmering 


treedt enige resolutievermindering op en een vervaging 
van de contouren, Dit is het gevolg van een te beperkt 
model. Het beeld is namelijk in zijn geheel beschreven 
m.b.v. een homogeen model, terwijl contouren en sprongen 
In een beeld kunnen worden gezien als gebieden waar het 
model verandert, Een mogelijke oplossing voor dit pro- 


bleem Is gegeven door (Biemond, Gerbrands, 1979). 


Vervolgens wordt onderzocht in welke mate zowel ver- 
vorming als ruis uitgefilterd kan worden. We doen dit 
eerst aan de hand van een eenvoudig door de computer ge 
genereerd J2xj2 testbeeld, voorstellend de letter "'G" 
(Figuur 12, linksboven), In de figuur rechtsboven is een 
lineaire bewegingsvervorming geïntroduceerd in horizon- 
tale richting over 6 beeldpunten (uitsmering over 6 beeld- 
punten), waarna additieve ruis wordt toegevoegd met 
SNR = 100 (linksonder). Tenslotte is in de figuur rechts- 
onder het filterrestultaat weergegeven. Hetzelfde expe- 
riment is herhaald voor het báx64 testbeeld van het 
meisjesgezicht (figuur 13), Duidelijk is te zien dat inhet 
gefilterde beeld de bewegingsvervorming nagenoeg is ver- 
dwenen en het beeld zijn oorspronkelijke scherpte terug- 


krijgt. 


Fig. 12: Filterresultaat 32x32 testbeelden, 
lineaire bewegingsenscherpte over 6 beeldpunten: 


SNR = 100, 


Fig. 13: Filterresultaat báx64 testbeelden, 
lineaire bewegingsonscherpte over 6 beeldpunten; 


SNR = 100. 


In figuur 14 wordt tenslotte het resultaat getoond van 
ons laatste experiment, Het 32x32 testbeeld (linksboven) 
wordt gedefocussserd m‚b.v, het nearest neighbour model 
(rechtsboven) en voorzien van additieve ruis met SNR = 
100 (linksonder), waarna rechtsonder het filterresultaat 


is afgebeeld, 


Fig. 14: Filterresultaat 32x32 testbeeld, 
defocussering volgens nearest neighbour model; 


aNR = HO, 


De auteur hecht eraan, de samenwerking met K.K. Eerland, 
Ir. J.J. Gerbrands en Prof.ir. IJ, Boxma op deze plaats 


te noemen. 
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Voordracht gehouden op 30 oktober 1980 op de THE tij- 
dens een gemeenschappelijke vergadering van het NERG 
(nr.292), de Sectie voor Telecommunicatietechniek KIvI, 


en de Benelux sectie IEEE. 
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SECTIE TELECOMMUNICATIE, Klvl 
NEDERLANDS ELEKTRONICA- EN RADIOGENOOTSCHAP 
(293ste werkvergadering) 

IEEE BENELUX SECTIE 


UITNODIGING 
voor de lezingendag over „Alternatieve methoden van energie-opwekking” op 
donderdag 20 november 1980 in de grote zaal van het gebouw van het Klvl 
aan de Prinsessegracht 23 te ‘s-Gravenhage. 


PROGRAMMA 
9.30 uur: Ontvangst en koffie. 


10.00 uur: DRS. K, JOON, (Energie-onderzoek Centrum Nederland, Petten): 
HET NATIONALE ONDERZOEKPROGRAMMA ZONNE-ENERGIE. 


10.30 uur: IR. C. DEN OUDEN, (Technisch Physische Dienst-TNO, Delft): 
THERMISCHE TOEPASSINGSMOGELIJKHEDEN VAN ZONNE-ENERGIE IN 
NEDERLAND. 


11.15 uur: Koffiepauze. 


11.30 uur: DR. IR. A. KIPPERMAN, (Technische Hogeschool, Eindhoven): 
TOEPASSING VAN ZONNECELLEN, ENIGE SYSTEEMBESCHOUWINGEN. 


12.15 uur: Lunchpauze, lopende lunch in de wandelgangen. 


13.45 uur: IR. G. G. PIEPERS, (Energie-onderzoek Centrum Nederland, Petten): 
KAN NEDERLAND VAN DE WIND LEVEN? 


14.15 uur: IR. W.J. DE ZEEUW, (Technische Hogeschool, Eindhoven): 
ALS NEDERLAND IN DE TOEKOMST VAN DE WIND ZOU LEVEN... 
ENIGE SYSTEEMBESCHOUWINGEN. 


15.00 uur: Theepauze 


15.15 uur: IR. J. OUDERLING, (Dr. Neherlaboratorium-PTT, Leidschendam): 
EFFECTEN VAN WINDENERGIE-SYSTEMEN OP TELECOMMUNICATIE- 
VERBINDINGEN, 


16.00 uur: Sluiting. 


Aanmelding kan geschieden vóór 14 november a.s. 

Tijdens de lunchpauze bestaat de mogelijkheid gebruik te maken van een lopende lunch in de 
wandelgangen van het Klvl-gebouw. 

De lunchkosten kunnen worden voldaan door het overmaken van f 12,50 op postgiro 576595 
tn.v. penningmeester Klvl Sectie T te Oosterhout, onder vermelding van „Lunch energiedag”. 
De betaling dient vóór 15 november te zijn ontvangen. 


Het Klvl-gebouw is als volgt te bereiken: 

Per trein naar Den Haag-C.S., to. het station iet een hertekamp „Koekamp", daarachter bevindt 
zich het Malieveld, links van Koekamp en Malieveld ligt de Prinsessegracht. (ca. 5 min. gaans 
van C.S.) 

Per auto komende vanaf de richting Gouda eindigt de autoweg op de Utrechtsebaan, aan het einde 
links (Zuid-Hollandlaan), aan het einde weer links. 

De Prinsessegracht ligt nu rechts van u. Parkeergelegenheid in de omgeving van het Malieveld of 


Fareeganage Di het Centraal. Statibe. Namens de samenwerkende verenigingen, 


IR. R. C. STRATO. 
Den Haag, oktober 1980. Telefoon 070 - 782640 


Twintigste plenaire vergadering van de Internationale Wetenschappelijke Radio 
Unie, Washington D.C., 10 = 19 Augustus 1981. 
F. Louis H.M. Stumpers 


Universiteit Utrecht en Nijmegen 


Philips Nat.Lab. (ret.}) 


Terwijl het Bestuur, het Executieve Comi- 
té en het Coördinerend Comité reeds van 7 tot 
9 Augustus bij elkaar komen voor diverse zake- 
lijke besprekingen, begint het wetenschappe- 
lijk gedeelte met de officiële opening op 
Maandag 10 Augustus 1981 in het Hyatt Regency 
Hotel in Waskington, waar alle vergaderingen 
worden gehouden. Diezelfde dag zullen de negen 
Commissies hun nieuwe vice-voorzitter kiezen, 
en hun plannen voor de komende drie jaar be- 
spreken. De Commissie voorzitters van de afge- 
lopen periode treden nu af, en worden opge 
volgd door hen vice-voorzitter. 

Om de eenheid van het wetenschappelijk 
onderzoek in U.R.S.I. te benadrukken worden 
drie conferenties gegeven over onderwerpen van 
algemene interesse, waarbij alle deelnemers 
verwacht worden, telkens 's morgens om 9 
uur : 12 Augustus professor Shearman: Radio 
science and vceanography, 17 Augustus profes- 
sor Gordon Solar power satellite and tele- 
communicatiors, en 19 Augustus professor Fried- 
man : High energy astronomy. Ook de uitrei- 
king van de van der Pol en Dellinger medailles 
en van de Appleton Award, waarbij de bekroon= 
de sprekers over hun werk zullen praten, be- 
looft van algemene interesse te zijn. 

Er worden vier tweedaagse open symposia 
georganiseerd. Op 1l en 12 Augustus Remote 
sensing, voorziter dr. J.0. Thomas, secretaris 
dr. D.L. Croon. De vier zittingen betreffen 
Remote sensing of the sea (including polar 
oceans), organisator dr. Valenzuela. Remote 
sensing of land (including snow/ice cover of 
land). Organisator dr. Gjessing. Remote sen- 
sing of the lower and middle atmosphere. Orga=- 
nisator dr. Croon. Image processing in remote 
sensing. Organisator dr. J. Thomas. Op 13 en 
14 Augustus Millimeter and submillimeter 
waves, voorzitter K.J. Button (M.I.T,.) met zit= 
tingen over radioastronomie (dr. Baars) over 
complexe dielectrische eigenschappen van vaste 
stoffen en vloeistoffen (dr. Afsar}) over me= 
ting van vermogens (Stauffer en Tait), over 
de vrije electron maser en het gyrotron (dr. 


Sprangle) en een algemene zitting (K.J. Button). 


Eveneens op 13 en 14 Augustus een open 


symposium over Interactie van electromagnetische 
golven en biologische systemen, voorzitter S.W. 
Rosenthal. Dit symposium wordt gecoördineerd 

met de jaarlijkse vergadering van de Bioelec=- 
tromagnetics Society van 10 tot 12 Augustus, zo= 
dat vijf dagen in onderlinge samenhang worden 
gevuld. De U,R.‚.S.I. dagen zullen zowel uitgeno- 
digde sprekers aan het woord zien, als aange- 
boden bijdragen. 

Het vierde open symposium betreft mathe- 
matische modellen van radio propagatie en het 
wordt georganiseerd door dr. J. Wait op 17 en 
18 Augustus. 

Het open symposium over Remote sensing 
trekt reeds zoveel belangstelling, dat wellicht 
meer tijd ervoor nodig zal zijn (of parallel- 
zittingen) . 

Gedeeltelijk parallel met de open symposia 
organiseren de negen Commissies van U.R.S.I. 
hun wetenschappelijke zittingen. Bij Commissie 
A (meetmethoden) zijn zittingen gewijd aan Me- 
tingen in electrische netwerken. Oppervlakte 
toleranties voor antennes met grote precisie 
(met B en J). Tijd en frequentie metingen. Sig= 
naal en ruis metingen (met E). De toepassing 
van de optische fiber bij metingen. Cryogene 
metingen. 

Commissie B (Electromagnetische golven) 
heeft zittingen over meetmethoden. Strooiing en 
buiging. Geleide golven. Methodes voor hoge 
frequenties. Niet lineaire effecten. Media van 
onregelmatige samenstelling. Overgangsverschijn= 
selen in velden. Inverse strooiing. Nieuwe me- 
thoden. Antennes. 

Commissie C (Communicatie-, informatie- en 
netwerktheorie) behandelt 
in ruis. Microelectronische netwerken en sys= 


Optimale codering 


temen (met D). Telecommunicatie en digitale sig=- 
naalverwerking. Hulpmiddelen voor zeer snelle 
verwerking van signalen (met D en J). Satelliet 
communicatie. Computer communicatie."Spread spec- 
trum" communicatie. Optische communicatie. Com- 
municatie en data verwerking voor ruimte onder- 
zoek. 

Commissie D (Electronics) Optische sys= 
temen en componenten. 

Commissie E (Ruis van aardse oorsprong en 
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door mensen gemaakte ruis. Effecten van niet=- 
=gaussische ruis op communicatie. 

Commissie F (Tropospheer) : Interferentie pro- 
blemen bij voortplanting en efficient gebruik 
van het radio spectrum. Strooimechanismen voor 
radio golven in de midden atmosfeer (met G). 
Optische propagatie. Propagatie tussen aarde 
en satelliet, systeem en radio meteorologische 
aspecten. Speciale onderwerpen. Rapporten over 
internationale symposia. 

Commissie G (Ionosfeer) : Equatoriale 
onregelmatigheden in de ionosfeer (met H). 
Aeronomische studies met behulp van digitale 
ionosferische sondes. Wijzigingen in de iono- 
sfeer. Invloed van de ionosfeer op radio sys- 
temen. Radio onderzoekingen van de ionosfeer 
op hoge breedten, eerste resultaten van 
EISCAT (met H). 

Commissie H (Plasmas) : Bepaling van 
plasma golf spectra op afstand. Injectie van 
magnetosferische golven. Plasma golf analyse 
met hulp van de computer. 

Commissie J (Radioastronomie) : Nieuwe 
ontwikkelingen in observatoria en laboratoria. 
Geschiedenis van de radio astronomie met spe- 
ciale aandacht voor de vijftigste verjaardag 
van Jansky's ontdekking, 

Zoals men ziet, werken verschillende com- 
missies samen bij de behandeling van onder= 
werpen. We hebben de zittingen dan alleen bij 
de eerst in het alfabet komende commissie 
vermeld. 

Men mag interessante voordrachten over 
een breed gebied verwachten. Het is te hopen, 
dat Nederland ondanks de ongunstige omstan- 
digheden goed vertegenwoordigd kan zijn. Even- 
tueel kan een gezamenlijke reis de kosten 
verminderen, 
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UIT HET NERG 


LEDENMUTATIES 


Voorgestelde leden 
Ir. J. Bruijn, Beekstraat 45, Apeldoorn 


Ir. J. van Cappellen, Vreelandseweg 58, 
Nederhorst den Berg 
Dr. K.S. Chung, Keversberg 4, Eindhoven 
Ir. H.T. Huele, IJsselstraat 17, Hengelo (0) 
Ir. M,.H.P. van Iersel, Kreitenmolenstraat 99, Udenhout 
Ir. A.C, Sonneveldt, A. van Scheltemaplein 62, Delft 
Ing. C.P.M. Willems, Immenhof 30, Teteringen 
Ir. 5S.C.J. Worm, Goorland 6, Bladel 


Nieuwe leden 

J.M.H. Hanselmen, Morgental 20, Krimpen a/d IJssel 
Ir.C.G.M. van 't Klooster, Vestdijkhove 68, Zoetermeer 
Ir. C. Noort, Anne Franklaan 189, Bussum 

Ir. J.P. Schuddemat, Roskam 28, Beuningen (Gld) 


Nieuwe adressen van leden 
R.W. Budding, Berlagelaan 139, Bilthoven 
Ir. A.A. Dogterom, Eikenlaan 11, Hilversum 
Ir. R. Dufour, Hagelkruisweg 75, Deurne 
Ir. P‚J. Haubrich, Alpenlandstrasse 1À, 
A2380 Perchtoldsdorf, Oostenrijk 


Ir. M.E. Kopijn, Dingenaarserf 6, Eemmes 
Ir. A. Mawira, Burg. Caan van Necklaan 259, 
Leidschendam 
Ir. A.H. Toet, Koetlaan 25, Delft 
Ir. J.A. Verhoef, Molendijk 34, Hamont Achel (Limburg), 
België 


Voorlopig progremma van werkvergaderingen 


9/4 Jaarvergadering.Patroonherkenning 

23/4 IC's T.H. Twente 

1/5 Digitale transmissie. Dr.Neherlaboratorium 

2/9 Elektronische filters 

29/9 Het Instrument. Thema nog niet bekend. RAI Amster- 
dam 

22/10 100 jaar telefonie RAI Amsterdam 

?/10 Digitale lokale netten 

2/11 Elektronica in lucht en ruimtevaart 


rr Telecommunicatie en energie 


VARIA 


Zoals in de algemene vergadering op |&-4-1980 is be- 
sproken, is de contributie m.i.v., 1981 vastgesteld op 
f 60,-- per jaar. 

Voor diegenen die de contributie voldoen via de 
periodieke overschrijving zal de wijziging collectief 
aan de girodienst worden doorgegeven. Eventuele bezwa- 
ren kunnen bij de penningmeester NERG, postbus 39, 
2260AA Leidschendam worden ingediend. 
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NEDERLANDS ELEKTRONICA- EN RADIOGENOOTSCHAP 
(294ste werkvergadering) 
SECTIE TELECOMMUNICATIETECHNIEK, Klvl 


IEEE BENELUX SECTIE 
mmm 


UITNODIGING 


voor de lezingendag op donderdag 29 januari 1981 in het Auditorium van de Technische 
Hogeschool te Eindhoven. 


Thema: SATELLIET - OMROEP 


PROGRAMMA 
9.30 uur: Ontvangst en koffie. 


10.00 uur: PROF. DR. IR. J. J. GELUK, (Wereldomroep, Hilversum): 
DE HEMELSE OMROEP. 


10.50 uur: DR. G. RADIG, (Messerschmitt-Bölkow-Blohm, München): 
DESCRIPTION OF THE FRANCO-GERMAN TELEVISION BROADCAST 
SATELLITE PROGRAM TV-SAT/TDF-1. 


11.15 uur: Koffiepauze. 


11.45 uur: DR. D. FASOLD en L. HEICHELE (Messerschmitt-Bölkow-Blohm. München): 
ANTENNA SYSTEM OF THE GERMAN BROADCASTING SATELLITE. 


12.30 uur: Lunchpauze. 


14.00 uur: E. KAMMER, (Philips, Eindhoven): 
SATELLIET-OMROEP EN KABEL-TELEVISIE IN NEDERLAND. 


14.45 uur: Theepauze. 


15.15 uur: IR. J. MENDRIK, (Dr. Neher Laboratorium, Leidschendam): 
CENTRALE ANTENNE-INRICHTINGEN: NU EN IN DE TOEKOMST. 


15.45 uur: Discussie. 


16.00 uur: Sluiting. 


Aanmelding kan geschieden vóór 22 januari 1981 door inzending van de aangehechte kaart, 
ingevuld en gefrankeerd met een postzegel van 45 cent. 

De mogelijkheid bestaat om op de TH de lunch te gebruiken. Daartoe kunnen bonnen worden 
besteld door een bedrag van f 12,50 over te maken op postgiro 2389759 tnv. J. Neessen te 
Woerden onder vermelding van „Lunch SAT-dag”. De betaling dient te zijn ontvangen vóór 
22 januari 1981. 

Bestelde bonnen kunnen voor aanvang van de lezingen worden afgehaald bij de organisatietafel. 
Op de dag zelf zijn geen bonnen meer verkrijgbaar. 


Namens de samenwerkende verenigingen, 
IR. J. NEESSEN. 
Telefoon overdag 070 - 755591 

Woerden, december 1980. Telefoon ‘s-avonds 03480 - 14539 
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Inter-noise 81. 6-8 oktober 1981 Rai-congress 
Amsterdam. Ronne Inter-No Ae 81: P.0.Box 85542; 


Melecon '81. 24-27 mei 1981 — Tel Aviv Hilton. Kontakt- 
adres P.O.B. 3473, Tel Aviv 63431, Israël. 


an aen conference on digital satellite 


Eckene, OEfices me Fifth TEE: rnationa 

on Digital Satellite Communie: don I | 
zionale delle Vomanteazto vie ee Vilta 
Piaggio; 16125 Genova, Italy. 
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Telecommunication Transmission into the digital era; 
17-20 maart Daal Weeen NE Ee konten: 
Engineers: Ee en A NO Oers De jhone 
01-2401871 (Ext.) 261 Telex 261176 


First Conference on Integrated Optics; 14-15 septer 
1981 IEE, Savoy Place, London WC2R OBL, UK Kontakt 
adres: identiek aankondiging hierboven. 


control and its Applications; A Ee 


nference on hot carriers in 
cor Heen 7-10 a eN Montpellier, France; 


Prof. J.P. Nougier, Université de Sciences 
et Techniques de Languedoc; Place Eugéne Bataillon; 
34060 Montpellier France. Telefoon (67) 54.71.84 


B Doxlander Strasse 72; 7500 Pen Telefon 
(0721) 59665020 


Audio evolution in the 80s; 17-20 march 1981; Congress 
centrum Hamburg Am Dammtor; Kontakt adres: Miss T.K.S. 
Bakker; he Polygram B.V.; P.0. Box 23; 3740AA Baarn= 
Holland 
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